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LECTURE 13. INHOMOGENEOUS EQUATIONS 

We discuss various techniques for solving inhomogeneous linear differential equations. 

Variation of parameters: the Lagrange procedure. Let us consider the linear second-order differ- 
ential operator 

(13.1) Ly = y" + p(t)y' + q(t)y 

with variable coefficients. 

If a nonvanishing solution of a homogeneous equation Ly = is known, then the corresponding 
inhomogeneous equation Ly = f can be solved, in general, by two integrations. It was discovered 
by Lagrange that if two linearly independent solutions of Ly = are known, then the inhomoge- 
neous equation Ly = f can be solved by a single integration. 

Let u and v be a pair of linearly independent solutions of Ly = 0, and form the expression 

(13.2) y = au + bv. 

If a and b are constant, this represents the general solution of Ly = 0. We will the inhomogeneous 
equation Ly = f by choosing a trial solution of this form, but with a and b functions of t, rather 
than constants. The method is called the method of variation of parameters. 
Let a and b are differentiable functions of t. By differentiation, 

y = (au + bv') + (au + b'v). 

We require 

(13.3) a'u + b'v = 

so that y' = au' + bv' . This simplifies the calculation of the second derivative, and 

y" = (au" + bv") + (au + bv). 

Therefore, 

Ly = y + py' + qy = aLu + bLv + au + bv = a'u' + b'v' . 

The second equality uses that Lu = Lv = 0. 

Solving Ly = f in the form in (13.2) then reduces to the linear system 

au + b'v = 0, 
a u + b v = j, 

in the unknown a' and b'. In the matrix form, 

(u v\ 
\u' v'J 

By Cramers rule, we solve the system, and 
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Here, the notation | • | stands for the determinant of the matrix. The denominator is the Wromskian 
W(u, v), so that we may write them as 

/ -fv y fu 



a 



W(u,vy W(u,v) 
Finally, by integration, we obtain the Lagange formula 

(13.4) y(t) = u(t) [ ~/ V dt + v(t) [ T J U dt. 



W(u,v) J W(u,v 

Lagrange's procedure extends to equations of order n and it represents an important advance 
in the theory of differential equations. 

A similar idea already appeared. For example, when studying the linear first-order differential 
equations, we replaced the homogeneous solution ce p by ve p , where v is a function. 

Example 13.1. Consider the Euler equation 

(13.5) x 2 y" - 2xy' + 2y = x 2 f(x), x > 0, 

where the prime denotes the differentiation in the x-variable. 
By the technique discussed in the previous lecture, we compute 

u = x, v = x 2 , W(u, v) = x 2 . 



For x > 0, we write (13.5) as 
Then Lagrange's formula (13.4) gives 



2 2 

y — y + = 

X x z 



y(x) = —x / f(x)dx + x 2 



X 2 

Exercise. If f(x) = x m , where m is a constant, in the above example, show that a particular 
solution of (13.5) is 

— xlogx if m = — 1, 

x 2 log x if m = 0, 



x m+2 

otherwise. 



y P (x) = 

k m(m + 1 

The general solution of (13.5) is y(t) = c\x + cix 2 + y p (x). 

The Green's function: initial value problems. As an important application of the formula (13.4) 
we can find an integral representation of the initial value problem for Ly = f, where L is given in 
(13.1). 

Let to be a point on the interval /. Integrating (13.4) from to to t, 

v(t) =<t) -^^dt + v(t) - ww dt 

f l u{t')v{t)-u{t)v{t') 
J t0 u(t')v'(t')-u'(t')v(t') J{ > 

This function satisfies the conditions 

y(to) = 0, y'{t ) = 0. 
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Indeed, y(t) = a(t)u{t) + b{t)v{t) and y'(t) = a(t)u'{t) + b{t)v'{t) where a{t) = ^^y^ dt ' and 
b(t) = mu{t \ t> 

[ ) wit') M - 

In summary, the function defined as 



/,„ * ' u (t'y(f) -«'(*')»(«') 

solves the initial value problem, 

Ly = /, y(to) = 0, y'(io) = 0. 
The function G(t, if) is called the Green's function. 

Example 13.2. We continue studying the Euler equation (13.5) satisfying the initial conditions 

y(xo) = 0, y'(xo) = for some xq > 0. 
The solution has an integral representation 

f x f(t) 
y(x) =x (x- t) J -^-dt. 

Jx t 

For example, if f(x) = x sin x then 

y(x) = x (x — t) sin tdt = x(x — xq) cos xq — x(sin x — sin xq) . 

JXQ 

Exercise. (The Green's function: boundary value problem) We consider the boundary value prob- 
lem 

y" +p(t)y' + g(t)y = f(t) on (t u t 2 ), y ( tl ) = y(t 2 ) = 0. 

If u and v are linearly independent solutions of the homogeneous equation y" + py' + qy = 0, then 
show that the solution of the boundary value problem is given by 

y(t)= G(t',t)f(t')dt', 
hi 



where G(t',t)= < 



W(t') ' 
u(t)v(t') 

if i < t' ^ t 2 . 



The method of annihilators. We introduce another method of finding a particular solution of 
linear inhomogeneous differential equation with constant coefficients. Let 

Ly = + piy (n - 1} + • • • + p n _iy' + p n y, 

where pj are real constants. We study the differential equation Ly = f, where / is a sum of 
functions of type 

fe xt , fe^smut, tV* cos i/t. 

Note that these functions arise as basis solutions of linear homogeneous differential equations 
with constant coefficients. We find a differential operator A satisfies Af = 0, then we reduce 
solving Ly = f to solving the homogeneous equation LAy = 0. Such an operator A is called an 
annihilator of /. 

We illustrate with an example. 
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Example 13.3. We consider the differential equation 

(13.6) y" - by' - 6y = te l . 

Let L = D 2 - 5D - 6 = (D - 2){D - 3). Then (13.6) is written as Ly = teK 

By the exponential shift law for D, we recognize that te l is a solution of the differential equation 
(D - l) 2 y = 0. In other words, (D - l) 2 is an annihilator of te l . Applying (D — l) 2 in (13.6), we 
obtain the homogeneous differential equation 

(D-2)(D-3)(D-l) 2 y = 0. 

It is easy to see that e*, te t , e 2t , e 3t form a basis of solutions of the above equation. Hence, we set 
a solution of (13.6) as 

y(t) = cie* + c 2 te l + c 3 e 2t + c 4 e 3 *, 

and determine the constants Cj. 

Since Le 2t = and Le 3t = 0, moreover, we may set 03 = 04 = 0. Hence, 

y (t) = cie* + c 2 te t . 

We compute 

Ly = (D 2 -5D- 6)(cie* + c 2 te t ) = (2ci - 3c 2 )e* + 2c 2 te t = te l 
to obtain c\ = 3/4 and c 2 = 1/2. Therefore, a particular solution of (13.6) is y(t) = 3/4e* + l/2te t . 
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